The Field- Tuned Superconductor-Insulator Transition with and without Current Bias 
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The magnetic-field-tuned superconductor-insulator transition has been studied in ultrathin Beryl- 
hum films quench-condensed near 20 K. In the zero-current limit, a finite-size scaling analysis yields 
the scaling exponent product vz = 1.35 ± 0.10 and a critical sheet resistance Rc of about 1.2Rq, 
with Rq — h/Ae^. However, in the presence of dc bias currents that are smaller than the zero-field 
critical currents, pz becomes 0.75 ± 0.10. This new set of exponents suggests that the field-tuned 
transitions with and without dc bias currents belong to different universality classes. 

PACS numbers; 73.50.-h, 74.76.-w, 74.25.Dw, 74.40. +k 



The superconductor-insulator (SI) transition in ul- 
trathin films is believed to be a continuous quantum 
phase transition occurring at T = as the quantum 
ground state of the system is tuned by varying disor- 
der, film thickness, magnetic field, or carrier concentra- 
tion through a critical value. Models of the SI transition 
can be grouped into two categories: ones that rely on 
fluctuations in the phase p|-p[ and others that rely on 
fluctuations in the amplitude |o] of the order parameter 
to drive the transition. For models in the first category. 
Cooper pairs are assumed to be present on the super- 
conducting side as well as on the insulating side of the 
transition. Based on this "dirty-boson" model. Fisher 
et al. have proposed a phase diagram for a two- 

dimensional system as a function of disorder and mag- 
netic field, with the transition being described in terms 
of interacting bosons (the Cooper pairs) in the presence 
of disorder. In the vicinity of the quantum critical point, 
the resistance is predicted to obey the scaling law HH]: 




Ri6,T) RJ{St), 



(1) 



where t{T) = T^^^'^^ , and 5 is the deviation of a tuning 
parameter from its critical value. For field-tuned transi- 
tions, 5 = \B — Bc\ with B and Be being the applied mag- 
netic field and the critical field, respectively. The criti- 
cal resistance Rc is predicted to have a universal value 
of Rq = h/Ae^ « 6.5 kfl/O. Scaling arguments set a 
lower bound on the correlation length exponent > 1 0] 
and give the value of the dynamical critical exponent z 
~ 1 H,^. Experimentally, critical exponents consistent 
with these predictions have been found in scaling analy- 
ses of the field-tuned transitions in InOa; and MoGe 
l^. However, the predicted universal critical resistance 
is not observed in these experiments. 

In this Letter, we report studies of the field-tuned SI 
transition in ultrathin Be films quench-condensed near 20 
K. We describe scaling analyses of the field-tuned tran- 
sition based on the " dirty boson" model. We have found 
that, for measurements performed in the zero-current 
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FIG. 1. Sheet resistance versus temperature measured on 
one marginally superconducting film at various magnetic field 
values. For curves from bottom to top, the field increased 
from 0.05 T to 2.25 T. The critical field was 0.66 T. In 
zero-field, ~ 2 K. 



limit, the scaling exponent product vz — 1.35 ± 0.10. 
This agrees well with results of the field-tuned transition 
in relatively thicker InO^; and MoGe ||^ films, but 
disagrees with a recent report 10| of lyz ~ 0.7 for the 
field-tuned transition in ultrathin Bi/Ge films. For our 
Be films having a robust superconducting phase, Rc ap- 
pears to be near 1.2Rq in the field-tuned transition. This 
agrees well with results from the Bi/Ge films [^,|l^. We 
have also carried out for the first time studies of the field- 
tuned transition in the presence of a dc bias current. The 
applied dc bias current should exert a force of j x B on 
the vortices in the direction perpendicular to both the ap- 
plied magnetic field B and the current density j. We have 
observed that the scaling exponent product vz becomes 
0.75 ± 0.10 in the presence of dc bias currents, suggest- 
ing that the field-tuned transitions with and without dc 
bias current belong to different universality classes. 

In Fig. 1, we show, with zero bias current, the temper- 
ature dependence of the sheet resistance, Ra, measured 
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at various field values for one of our Be films with the 
magnetic field applied perpendicular to the plane of the 
film. These Be films were quench-condensed onto bare 
glass substrates which were held near 20 K during evap- 
orations, under UHV conditions inside a dilution refriger- 
ator. This in situ progressive evaporation setup allowed 
for systematic studies of the SI transition as film thick- 
ness was varied. We deposited each set of films in fine 
steps. We carefully monitored Ru during each evapo- 
ration step until a desirable value of Ru was reached. 
The films were very close to 10 A in thickness, however, 
our quartz thickness monitor was not sensitive enough to 
pick up the small thickness increments after each step. 
Film resistance was measured in a standard four-terminal 
geometry using a PAR-124A lock-in amplifier operating 
at 27 Hz. The ac probe current was fixed at 1 nA. At 
the finite measuring temperatures in the vicinity of the 
field-tuned transition, the J- F characteristics in this low- 
current regime were linear. More details regarding these 
Be films have been published elsewhere Quench- 
condensed Be films are chosen for this study because such 
films were found to be nearly amorphous [ pT[ . These Be 
films undergo a transition from insulating to supercon- 
ducting when the normal state sheet resistance, Rn^ 
is reduced below ~ 10 kJl/D with increasing film thick- 
ness. We have now studied the field-tuned SI transition 
in several films of Rn between 5.6 and 12 kO/D. The su- 
perconducting transition temperature, Tc, of these films 
varied between 0.5 to 4 K in zero-field. 

The film in Fig. 1, with Rn = 10.7 VVl/U at 15 K, is 
considered marginally superconducting, as we will discuss 
later. With increasing field, corresponding to curves from 
the bottom to the top in Fig. 1, this film was driven from 
superconducting to insulating, with a rather fiat Ra vs. 
T curve at a critical field of Be = 0.66 T. The main part of 
Fig. 2(a) shows Ru vs. B at various temperatures for the 
same film as shown in Fig. 1. In the vicinity of a quan- 
tum critical point, the resistance of a two-dimensional 
system is predicted to obey the scaling law in Eq. (1). 
Determining the critical exponents involves plotting the 
Ru vs. B data at various temperatures according to the 
scaling law. The good crossing point, over one decade in 
temperature, in the Ru vs. B plot in the main part of 
Fig. 2(a) identifies B^ = 0.66 T and R^ = 4.4 VVt/U for 
this film. We have used two methods to determine the 
scaling exponent vz. First, we can find vz by evaluating 
the derivative of the Ru vs. B curve at the critical field. 
In this case, we have the following scaling equation: 

— |B.cxi?,r-i/-/'(o). (2) 

A plot of vs. on a log-log scale, shown in 

the inset to Fig. 2(a), yields a straight line with a slope 
equal to l/vz, from which we determine vz — 1.36 ± 
0.10. Alternatively, we can plot R/Rc vs. \B — Bc\t 
and treat t(T) as an unknown variable. The values of 
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FIG. 2. Main figures show sheet resistance vs. field at a 
number of temperatures as labeled in the figures, with hias ~ 
in (a) and luas ~ 2.5 /xA in (b). The crossing points iden- 
tify the critical fields and critical resistances. Insets show the 
power-law relations between f§|sc and for the corre- 
sponding Ibias values, which determine the values of uz. 

t{T) at various temperatures are determined by obtaining 
the best collapse of the data. Following this procedure 
presented in recent papers ||9|,|l0|] by N. Markovic et al., 
we obtain the temperature dependence of t{T) which we 
plot on a log-log scale in the inset to Fig. 3(a). The 
collapse of the data is shown in the main part of Fig. 3(a) 
in a R/Rc vs. \B — Bc\t plot. Figure 3(a) shows good 
collapse of the data over three orders of magnitude m\B— 
Bc\t and two orders of magnitude in R/Rc- The straight 
line in the inset to Fig. 3(a) shows a power-law fit, as 
expected by the scaling function, Eq. (1). The slope of 
the line in the inset gives -l/vz, from which we find vz — 
1.34 ± 0.10. The exponents obtained from the above two 
methods agree with each other, showing the consistency 
of the scaling analysis. The scaling exponents obtained 
in our Be films in the zero-bias limit appear to be in very 
good agreement with the predictions of theories based 
on the "dirty boson" model Our result is also in 

very good agreement with renormalization group theories 
[|l2|,^ and Monte Carlo |g simulations. 

Another important prediction of the bosonic model 
is the universal critical resistance Rc = Rq = h/Ae^ 
K, 6.5 kfJ/D at the quantum critical point. This re- 
mains a controversial issue since only in the Bi/Ge sys- 



2 



100 




B-BJt 

FIG. 3. Main figures show the scahng plots of R/Rc vs. 
\B — Bc\t, with Ibias = in (a) and hias = 2.5 fiA in (b). 
Insets show the power-law relations between the parameter t 
and the temperature, which determine the values of i/z. 



tern I^JT^ has a critical sheet resistance close to the pre- 
dicted value of Rq been observed. In fact, careful in- 
vestigations in the Bi/Ge system have revealed that Rc 
is l.li?Q to 1.2Rq and decreases in a narrow range as 
Rn is reduced with increasing film thickness. It was sug- 
gested that the slightly larger Rc than Rq was con- 
sistent with calculations in two-dimensions for a bosonic 
model including Coulomb interactions |l4| , which pre- 
dicted Rc ^ IARq, as well as Monte Carlo simulations in 
the (2+l)-diniensional XY model without disorder [p^ , 
which found Rc = 1.2Rq. The small variation of Rc with 
Rn could be explained pO(| as a finite-temperature effect, 
since strictly speaking the critical resistance is predicted 
to be universal only at T = 0. In other systems, such as 
MoGc films, Rc varies in a much wider range among 
films of varying i?jv- In order to explain the low Rc of 
their MoGe films, Yazdani and Kapitulnik proposed a 
two-channel model in which the conductance due to nor- 
mal electrons add to the conductance due to the Cooper 
pairs. If the unpaired electrons are localized, this model 
collapses to a single channel bosonic model with Rc — 
Rq. Otherwise, unpaired electrons add to the conduc- 
tion at the transition. It has been noted that the 
resulting film resistance can be either larger or smaller 
than Rq since the magnetoresistance contribution p^ ] 
from the unpaired electrons can be either positive, if the 
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FIG. 4. Figures show (a) the critical field. Be, and (b) the 
critical resistance, Rc, as functions of the normal-state sheet 
resistance, Rn, measured at 15 K. 



spin-orbit coupling is strong, or negative, if the spin-orbit 
coupling is very weak. 

Our Be films showed a zero-field SI transition as Rn 
was reduced below '^10 kSl/D with increasing film thick- 
ness 0. Films of Rn near 11 kf2/n are considered 
marginally superconducting since they have near-zero 
critical field values. Further increasing film thickness, 
the critical field. Be, increases with decreasing Rn, as 
shown in Fig. 4(a). In Fig. 4(b), we plot how Rc in the 
field-tuned transition varies with Rn ■ For marginally su- 
perconducting films of Rn between 9 to 12 kfi/D, Rc 
was significantly smaller than Rq. Nevertheless, we see 
clearly by comparing Fig. 4(a) with Fig. 4(b) that, as 
Rn was reduced by increasing film thickness, films with 
robust critical field values have a Rc of 7.1 ^ 8.0 kil/D 
in the field-tuned SI transition, which is about 1.2Rq. 
This result is in excellent agreement with results from the 
Bi/Ge system |10|] . However, our results disagree with 
the suggestion Q that the magnetoresistance of the un- 
paired electrons caused the discrepancy between Rc and 
Rq . Since Be has the weakest spin-orbit coupling among 
metals, the magnetoresistance should be negative, lead- 
ing to Rc < Rq following the two-channel model. We 
note that we have shown in Fig. 1, Fig. 2(a), and Fig. 
3(a) data measured on a marginally superconducting film 
of Rn = 10.7 kfi/D at 15 K. Although the critical resis- 
tance of this film, Rc ~ 4.4 kfl/D, is significantly smaller 
than Rq, we have found that, for all the marginally and 
robust superconducting films of Rn ranging from 5.6 to 
12 kfi/D, the critical exponents are the same with i^z = 
1.35 ± 0.10 in the zero-current limit. 

Below, we describe results of the field-tuned SI tran- 
sition in the presence of a dc bias current, hias- For 
such studies, Ibias was varied between 125 nA and 2.5 
fiA and kept below the zero-field critical current (~ 15 
fiA) . At each fixed Ibias , we used a magnetic field to tune 
a film from superconducting to insulating. We believe 
that joule heating was insignificant in our experiments 
based on the following arguments. First, we have per- 
formed extensive /- V measurements over the entire tem- 
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perature and magnetic field range of our experiments, 
with the time span for the /- V sweeps ranging from ten 
minutes to one hour. The /- V curves were completely 
reproducible without any observable hysteresis and inde- 
pendent of the sweep rate. Secondly, we can estimate 
the temperature increase, AT, on the Be films due to 
joule heating. The thin Be films were deposited on glass 
microslide substrates of thickness 0.23 mm, which was 
attached to a copper sample holder by a very thin layer 
of grease. The glass substrate was the dominant source of 
heat resistance, with a thermal conductivity of ^ 0.0003 
W/Km at 100 mK jl7j. For a typical film square of size 
3x3 mm^ and resistance 10 kfi, joule heating for Ibias = 
125 nA is about 0.15 nW, resulting in AT - 0.012 mK 
at 100 mK. For I^as = 2.5 ^J,A, AT - 5.0 mK at 100 
mK. We note that the heating power at 2.5 fiA was 400 
times larger than the heating power at 125 nA. If joule 
heating were significant, the data obtained with Ibias = 
2.5 fiA should show a flattening of the data in the low 
temperature region when compared to the data obtained 
with Ibias = 125 nA. The fact that the scaling results at 
125 nA and 2.5 /iA agree well suggests that heating was 
insignificant. This also argues against the existence of 
significant electron heating decoupled from the lattice. 

We plot in Fig. 2(b) and Fig. 3(b), for h^as = 2.5 ^J,A, 
the results of scaling analyses based on the two previ- 
ously described methods. In Table I, we list the parame- 
ters from scaling analyses of the field-tuned transition at 
various hias- Results from the data collapsing method 
are presented for Ibias values of 250 nA and 2.5 /iA, for 
which the amount of data taken was adequate for such 
analyses. It appears that vz ^ 0.75 ± 0.10, showing no 
systemic change with Ibias in the range we have studied. 
Nevertheless, it is significantly smaller than the vz found 
in the zero-current limit. We can only speculate that the 
bias current could lead to a symmetry-breaking, resulting 
in different critical behavior at the transition. We note 
that experiments need to be carried out in other systems 
in order to find out whether this result is universal. In ad- 
dition, experiments at low Ibias values should be carried 
out to determine whether vz changes abruptly or grad- 
ually as Ibias is increased from zero. Such experiments 
can probe the threshold Ibias for the change in the scaling 
exponents, and have the potential of revealing whether 
the new scaling exponents are produced by certain non- 
linear effects in the vortices under a bias current. On 
the other hand, we need to discuss other possible origins 
that might have led to the new scaling exponents. For 
example, transport in marginally superconducting films, 
or, in the presence of a dc bias, could be dominated by 
narrow superconducting filaments, leading to a change 
in the dimensionality of the system. While it is difficult 
to determine to what extent film inhomogeneity affects 
our experiments, we would like to comment on the crit- 
ical currents measured in our films. For the film shown 
in Fig. 2(b) and Fig. 3(b), the critical current in zero 



TABLE I. Bias currents and derived parameters 



Ibias (nA) 


Be (T) 


Rc {kn/a) 


Rn (kn/D) 


vz * 


4. 

uz ' 


125 


0.33 


15.5 


9.36 


0.81 




250 


0.35 


16.1 


9.36 


0.76 


0.73 


1000 


0.36 


15.5 


9.36 


0.73 




2500 


0.38 


15.0 


9.36 


0.77 


0.75 



*i'z obtained by the fjlsc method 



^I'z obtained by the data collapse method 

field was about 1.5x10^^ A. This film was about 10 A in 
thickness and 3 mm in width. Thus, the critical current 
density would be about 5x10^ A/m^ if the current were 
uniformly distributed in the bulk of the film. On the 
other hand, if we assume that the current runs through a 
few filaments of 100 A in total width, the critical current 
density would be 1.5x lO^^A/m^. The critical currents of 
amorphous Be films have been measured by other groups, 
for example by Okamoto et al. jl^. Their Be films had 
a resistivity of 3.6x10"^ Vim. This value gives rise to 
Rn ~ 3.6 kil/D for a film of thickness same as ours (« 
10 A), meaning that their films were roughly 3 times less 
resistive than ours. They found that their films had crit- 
ical current density of about 1x10^ A/rs?. Compared to 
their critical current density, the assumption in our films 
of conduction dominated by a few narrow filaments leads 
to an unreasonably large critical current density. In ad- 
dition, we can compare our marginally superconducting 
films with robust superconducting films which are un- 
likely to be dominated by narrow filaments. As we have 
discussed, with Ibias = 0, these two types of films showed 
identical scaling exponents vz. It is possible that film in- 
homogeneity has lead to a depressed Rc in marginally 
superconducting films, but has not changed the effective 
dimensionality of the system. Therefore, the change in 
the scaling exponents with a dc bias is unlikely due to a 
dimensional cross-over as a result of filamentary super- 
conductivity. 

In conclusion, the scaling exponents found in our Be 
films in the zero-biased, field-tuned SI transition agree 
very well with results from InO^, and MoGc films, 
but disagrees with recent results |H,|l^ from quench- 
condensed Bi/Ge films. For Be films having a robust 
superconducting phase, the critical sheet resistance in 
the field-tuned SI transition was about \.2Rq. Our field- 
tuned transition in the presence of a dc current reveals 
a set of new scaling exponents, suggesting that the field- 
tuned transitions with and without dc bias belong to dif- 
ferent universality classes. We gratefully acknowledge S. 
Teitel and Y. Shapir for numerous useful discussions. 
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